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Abstract 



We summarize recent results obtained in (TJ on the running m c ^(m c ^) in the MS scheme and fD (s) ,B (s) using QCD spectral sum 
rules (QSSR) known to N2LO PT series, including all dimension-six NP condensate contributions in the full QCD theory, an 
estimate of the N3LO terms based on the geomteric growth of the PT series and using the most recent values of the QCD input 
parameters given in Table [T] The study of the effects of the subtraction scale [i on "different QSSR data" and the use (for the first 
CN time) of the Renormalization Group Invariant (RGI) scale independent quark masses in the analysis are emphasized. The estimates 
^ [rigourous model-independent upper bounds] reported in Table |2j f D /f„ = 1.56(5)[< 1.68(1)], /^/A = 1.58(5)[< 1.80(3)] and 
O f D /f K = 1.58(4)[< 1.63(1)],/^//* = 1.50(3)[< 1.61(3.5)], improve previous QSSR estimates. The remarkable agreements with 
some of the experimental data on f D D and with lattice simulations within dynamical quarks confirm both the success of the QSSR 
semi-approximate approach based on the OPE in terms of the quark and gluon condensates and of the Minimal Duality Ansatz 
jjy (MDA) for parametrizing the hadronic spectral function which we have tested from the complete data of the J I if/ and T systems. 
The running quark masses m c (m c ) = 1286(66) MeV and m^(m^) = 4236(69) MeV from M D ^ B are in good agreement though less 



CO 



accurate than the ones from recent J j 'if/ and T sum rules. 
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1. Introduction and a short historical review 

The (pseudo)scalar meson decay constants fp are of prime in- 
terests for understanding the realizations of chiral symmetry in 
QCD. In addition to the well-known values of f n and fx which 
control the light flavour chiral symmetries, it is also desirable to 
extract the ones of the heavy-light charm and bottom quark sys- 
tems with high-accuracy. These decay constants are normalized 
through the matrix element: 



(0\J* Q (x)\P) 



fpMj 



J p qQ (x) = (m q + M Q )q(iy 5 )Q , (1) 



X 



where: J?q(x) is the local heavy-light pseudoscalar current; 
q = d,s; Q = c, b\ P = D( S ),B( S ) and where fp is related to 
the lep tonic width T(P + — > / + v/). Since the original works of 
NSV2Z Q based on QCD spectral sum rules (QSSR) E)[] a 
large amount of QSSR works have been devoted to give bounds 
on fp tSHlQl and to estimate their values (TTUT3ll . More re- 
cent works including N2LO PT corrections have been derived 
later on in full QCD lfT4lfT6l and in HQET HZ). Different ear- 
lier papers (HE) have been scrutinized in 0[TT]], where Nari- 
son found that the apparent discrepancies between the different 
QSSR results can be solved if one applies carefully the stabil- 
ity criteria (also called sum rule window) of the results versus 
the external QSSR variables and continuum threshold t c . In 
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this way, and for given values of m c ^, he obtained the values: 
f D * (1.31 ± 0.12)/^ , f B (1.6 ± 0.1)/* , which are indepen- 
dent of the forms of the sum rules used. However, the result 
has been quite surprising as it indicates a large violation of the 
heavy quark symmetry scaling predictions (T8), where 1/Mq 
corrections have been estimated numerically in [12]. This "un- 
expected result" has been confirmed few years later by lattice 
calculations lfT9l . In this talk, we present improved estimates 
and bounds on f P and on m c ,b(^c,b) obtained recently in Q 
in full QCD theory, where the most recent values of the (non- 
)perturbative QCD parameters given in Table [TJhave been used. 
The expressions of NLO PT in OHO), of N2LO PT in ED, of 
the non-perturbative in (71 [10 ] and the light quark mass correc- 
tions in (91 OH |22l have been used. The N3LO PT contribu- 
tions have been estimated by assuming the geometric growth 
of the series lf23l which is dual to the effect of a l/q 2 term 
(24l [23. Like in fT3lL the previous PT results obtained in the 
on-shell scheme is translated to the one in the MS -scheme by 
using the relation between the on-shell and MS mass known at 
present to N3LO (HE1. The Renormalization Group Invariant 
(RGI) s, c, b quark masses introduced by (26l which are scale 
and (massless) scheme independent have been also used for the 
first time, while a careful study of the effect of the substraction 
scale on each "QSSR data point" has been performed. 

2. QCD spectral sum rules (QSSR) 

Forms of the sum rules 

We shall be concerned with the two-point correlator : 

<e^ 2 > = ' f d * x e iqx (0\Tjr Q (x)jr Q (0?\0) , (2) 
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where JqQ(x) is the local current defined in Eq. |T]). The as- 
sociated Laplace sum rules (LSR) X^q(t) and its ratio ^q(t) 
read OR 



- f 

J(m 



dt e 



1 



Imi// - Q (t,fi) , 



(3) 



where [i is the subtraction point which appears in the approxi- 
mate QCD series when radiative corrections are included. The 
ratio of sum rules ^g(r,yu) is useful, as it is equal to the res- 
onance mass squared, in the Minimal Duality Ansatz (MDA) 
parametrization of the spectral function: 
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Im^ Q (0 ^ fpM 4 P S(t - M 2 P ) + "QCD cont."6>(7 - t c \ (4) 



where f P is the decay constant defined in Eq. ([I]) and the higher 
states contributions are smeared by the "QCD continuum" com- 
ing from the discontinuity of the QCD diagrams and starting 
from a constant threshold t c . We shall also use for the 5-meson, 
the moments obtained after deriving n + 1 -times the two-point 
function and evaluated at Q 2 = (MSR) 13, where an expan- 
sion in terms of the on-shell mass M\, can be used. They read: 



A<>> 



Jf 

Jim 



(m q +M b ) 

■tc 



2 V 



dt 1, 

n+2 



(5) 



that the MDA, with a value of ^T c around the value of the 1st 
radial excitation mass, describes quite well the complete data in 
the region of r and n where the corresponding sum rules present 
r (in units of (GeV -2 ) or n stability 1281 : 

* (1.3 - 1.4) , r T * (0.2 - 0.4) , n r * (5 ~ 7), (6) 

as we shall see later on. Moreover, MDA has been also used 
in |29ll (called Minimal Hadronic Ansatz in this paper) in the 
context of large N c QCD, where it provides a very good ap- 
proximation to the observables one compute. 

a) 
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Figure 1 : a) r-behaviour of the ratio of £ff/£f b ml for y^T = M T(25) -0. 15 GeV. The red 
dashed curve corresponds to the strict equality for all values of r.; b) the same as a) but for 



Test of the Minimal Duality Ansatz (MDA) from J y 'if/ and T 
The MDA presented in Eq. Q, when applied to the p-meson 
reproduces within 15% accuracy the ratio measured from 
the total cross-section e + e~ — > I = 1 hadrons data (Fig. 5.6 of 
1 5 ]), while in the case of charmonium, M 2 from 9^ has been 
compared with the one from complete data where a remark- 
able agreement for higher n > 4 values (Fig. 9.1 of Q) has 
been found. Tests of MDA from the J I if/ and T systems have 
been also done in Q. Taking ^[T c My^s) - 0.15 GeV, we 
show (for instance) the ratio between and £!^ al in Fig. 

S) for LSR and M^ exp and M ( £ dual for MSR in Fig. jit) 
the T systems indicating that for heavy quark systems me 
role of the QCD continuum will be smaller than in the case 
of light quarks and the exponential weight or high number of 
derivatives suppresses efficiently the QCD continuum contri- 
bution but enhances the one of the lowest ground state in the 
spectral integral. We have used the simplest QCD continuum 
expression for massless quarks to order a 3 s from the threshold 
t c EUEl QCD cont. = 1 + as + 1.5as 2 - 12.07 as 3 . One can see 



2 Radiative corrections to the exponential sum rules have been first derived 
in [27), where it has been noticed that the PT series has the property of an 
Inverse Laplace transform. 

3 We have checked that the spectral function including complete mass cor- 
rections give the same results. 
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Figure 2: a) r-behaviour of ^(r) normalized to the ground state energy Eq for the 
harmonic oscillator. 2 and 4 indicate the number of terms in the approximate series. 

Optimal results from stability criteria 

Using the theoretical expressions of £!^ Q or M^ th , and 

parametrizing its experimental side or M^ exp by the 

MDA in Eq. ([4]), one can extract the decay constant f P and 
the RGI quark mass yHq. In principle the equality £f^ Q = £ e -* p 
should be satisfied for any values of the external (unphysical) 
set of variables (r, t c ), if one knows exactly and -Cj^. Here, 
we make an analogy with the harmonic oscillator, where the ra- 
tio of moments H^q is a function of the imaginary time variable 
r and where one knows the exact and approximate results. One 
can find l30l that the exact energy Eq of the ground state can 
be approached from above by the approximate series (see Fig. 
[2]). At the minimum or inflexion point (stability) of the curves, 
one has a ground state dominance. For small time (large Q 2 ), 
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all level contributes, while for large time (small Q 2 ) the series 
breakdown. We shall apply this stability criterion inspired from 
quantum mechanics in our analysis. In principle, the continuum 
threshold in Eq. ^ is a free parameter, though one expects 
its value to be around the mass of the 1st radial excitation be- 
cause the QCD spectral function is supposed to smear all the 
higher state contributions in the spectral integral as explicitly 
shown previously in Section [2] In order to avoid the model- 
dependence on the results, Refs. l4l l5l fTTUT4l have considered 
the conservative range of ^-values where one starts to have r- 
or ft- stability until which one reaches a ^-stability where the 
contribution of the lowest ground state to the spectral integral 
completely dominates. For the D and B mesons, this range is: 
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Figure 3: a) r-behaviour of fo from Xj c for different values of t c , for a given value of 
the subtraction point /i = r _1/2 and for m c - 1467 MeV; b) the same as a) but for M D from 



3. The QCD input parameters 



We shall use the QCD parameters (with obvious nota- 
tions) and their values given in Table [TJ where p ^ 2 
indicates the deviation from the four-quark vacuum satura- 
tion. From the running m q ,Q quark parameters, one can 
deduce the corresponding RGI quantities m qt Q and jl q l26l 
known to order a] EHUEU: ™> q ,Q(f) = m q ,Q(-Pia s )~ 11131 x 
C{a s \ (qq)(r) = (-fra^ /C(a s ), ' (qgoGqXr) = 

-Mfrli-fra^/Cias) , where A = -(1/2)(11 - 2n f /3) 
is the first coefficient of the J3 function for rif flavours; a s = 
a s {r) 1 7i. The QCD correction factor C(a s ) is 1 + lA155a s + 
1.5008^ + for nf = 5 flavours and shows a good con- 
vergence. Therefore, the RGI quantities to order a 2 s (heavy 
quarks) and to order a s (light quarks), in units of MeV are: 
fh c = 1467(14), m b = 7292(14), m, = 128(7), ft q = 251(6) . 
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Table 1: QCD input parameters. 



Parameters 



Values 



a s (M T ) 
m s {2) 
m c (m c ) 
m b (m b ) 



\{uu + dd) ll \2) 

(ss)l(dd) 

Ml 

(a s G 2 ) 
(g 3 G 3 ) 
p(qq) 2 



Ref. 



0.325(8) 
96.1(4.8) MeV 
1261(12) MeV 
4177(11) MeV 
-(275.7 + 6.6) MeV 
0.74(3) 

(0.8 + 0.2) GeV 2 
(7 + 1) x 10"2 GeV 4 
(8.2±1.0)GeV 2 x(of,G 2 > (28) 
(4.5 + 0.3) x 10- 4 GeV 6 l32l [371140 



E2HS1 

average t35l 
average lf28l 
average j|28l 

EES 

1 4 35][36l 
(37H39) 

|28l[30l[32ll40l[4n 
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Figure 4: a) n behaviour of f B from MSR for different values of t c , for m b =7292 MeV 
and at the subtraction point \i = 4 GeV; b) the same as a) but for M B . 



4. QSSR analysis and results 

The QCD expressions of the sum rules are given in Q to 
order a 2 s . We illustrate the analysis for the LSR in Fig. [3] for 
the D-meson and for the MSR in Figj4] for the 5-meson. The 
optimal values of f D(s) and f B(s) are obtained in the range of r and 
t c values given in Eqs. ^ and ([7]). We study the dependence of 
the results on the values of the subtraction scale [i in Fig. ([5]). 
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Figure 5: a) Values of fo from LSR at different values of the subtraction point fi and for 
ra c =1467 MeV: b Values of fn from LSR (red triangle) and from MSR (blue open circle) 
at different values of the subtraction point fi and for rhb=1292 MeV. 



Id - Ib have been estimated in |[T2l . f D and f Ds agree within 
the errors with the data compiled in (47H48), while the upper 
bound on f Ds can already exclude some existing data and theo- 
retical estimates (see e.g. (49 ] for some attempts to extract \V c d\ 
and \ V CS \)- Our results are compatible (in values and precisions) 
with lattice simulations including dynamical quarks I43U461 , 
which are not surprising as both methods study the same two- 
point correlator though evaluated in two different space-times 
and use the 1st principles of QCD (here is the OPE in terms 
of the quarks and gluon condensates which semi- approximate 
non-perturbative confinement). These agreements also confirm 
the accuracy of the MDA for describing the spectral function 
in the absence of complete data, which has been tested earlier 
IU |3 and in this paper from the charmonium and bottomium 
systems and in [ 29 ] in the large N c limit of QCD . 
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5. Summary and conclusions 

Table 2: Results from the open charm and beauty systems in units of MeV and 
comparison with experimental data and lattice simulations using rif = 2 143 11441 
and nf = 3 |45 46] dynamical quarks, fp are normalized as f n = 130.4 MeV 
mdf K = 156.1(9) MeV [47]. 



Charm 


Bottom 


Ref. 
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This work 
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f B 
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206(7) = 1.58(5)/, 


This work 


< 218.4(1.4) = 1.68(1)/, 
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fD, 


f Bs 




246(6) = 1.59(5)/* 


234(5) = 1.51(4)/, 


This work 


< 253.7(1.5) = 1.61(1)/* 
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260(5.4) 




Data 14711481 


248(6) 


232(12) 


ETMC(43) 
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26 Wj^ave, re,-e,xtracte,H fW£tffi)ay constants f D Jtyhk \ ^ a nd 
the running quark masses m c ^(Jn c ^) using QCD spectral sum 
rules (QSSR), the recent values of the QCD (non-)perturbative 
parameters given in Table [T] and (for the first time) the scale in- 
dependent Renormalization Group Invariant (RGI) heavy quark 
masses in the analysis after translating the on-shell PT expres- 
sion of the spectral function to the MS scheme. We have no- 
ticed that fp are very sensitive to mq and decreases when yuq 
increases. We have taken the conservative range of t c in Eq. 
([7]) covering the beginning of the r- or ft- stability until the te- 
stability [Figs. ^ and Q]. We have carefully studied the ef- 
fects of the subtraction scale \i on the "QSSR data" [Fig. (|5J]. 
Our final results in Table [2] agree and improve existing QSSR 
results in the literature. Large mass corrections responsible of 
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